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Abstract 

We study the coexistence of phases in a two-species model whose free energy is given by the scaling 
limit of a system with long range interactions (Kac potentials) which are attractive between particles 
of the same species and repulsive between different species. 
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1. Introduction. 

The rigorous ab initio derivation of the full phase diagram of even the simplest realistic model of a 
physical system, is still far beyond our mathematical grasp. This is so even when we restrict ourselves 
to simple classical one component system, e.g., point particles with Lennard-Jones type pair interactions. 
Only the low density-high temperature properties of such a system can be fully derived from statistical 
mechanics: there one can prove analytic behavior of the free energy as a function of the fugacity (or density) 
and temperature. Beyond that we have almost no means to prove the commonly observed facts of the phase 
diagram such as the Gibbs phase rule, the smoothness of the boundaries between the phases, etc. [17] 

To overcome this deficiency one can consider even more simplified models such as lattice systems. There, 
thanks to the Pirogov-Sinai [15], [22] theory, some exactly solvable models, inequalities, etc., one has much 
better control over the low temperature region of the phase diagram [17], [20], [7], [19]. Alternatively, one 
can consider, following the pioneering work of van der Waals [21], [18], [16], [5], mean-field-theories (MFT) 
yielding approximate equations for state or free energies which indeed exhibit most of the qualitative and 
many quantitative features of real world phase diagrams. These results were initially obtained in a heuristic 
manner and had to be supplemented by additional rules, e.g. the Gibbs double tangent or Maxwell equal 
area rule, to make them thermodynamically consistent. Subsequently, following the work of van Kampen 
[11], Kac, Uhlenbeck and Hemmer [9], Lebowitz and Penrose [12] were able to derive generalized mean field 
models in a more consistent mathematical way by considering systems in which part of the interactions 
were explicitly modeled by Kac potentials. Kac potentials are simply potentials of the form ^U^x), which 
contain a range parameter 7 -1 [9], [12]. For fixed 7 > these interactions are essentially finite ranged but 
in the limit 7 — > they become mean field like [12]. Depending on whether 7 -1 is small or comparable to 
the macroscopic scale (but always very large compared to the interparticle distances) one gets either the 
usual MFT (including the supplemental rules) or a macroscopic continuum theory (MCT) from which one 
can also obtain the surface tension associated with phase transitions caused by these Kac potentials [8] , [1] . 
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More recently it has been possible to prove phase transitions close to the mean field ones for large finite 7 1 
[13]. 

The structure of the coexisting phases in the MCT for a one component system has been extensively 
investigated in recent years by Presutti et al [6] and references therein. In multicomponent systems the 
MCT is rather unexplored. In the present work we study binary mixtures where the phase diagram has new 
qualitative features, resulting from the possibility of segregation of the system into regions rich in one of the 
species, e.g. oil and water. Such systems were already considered by van der Waals [21] and by Korteweg 
[10] using MFT expressions for the free energy. MFT however neglects the spatial structures of the phase 
domains so the starting point of our analysis is a MCT expression for the free energy J 7 of a binary mixture 
in a macroscopic domain A, which we shall take for simplicity to be a d-dimensional torus of side length L 
and volume |A| = L d . 

The free energy functional that we consider satisfies certain rearrangement inequalities, and we study the 
consequences of these for the phase diagram. We will later focus on special cases of physical interest, such as 
a mixture of van der Waals gases with a repulsive interaction between the species. However, to make clear 
the role of the rearrangements, we begin in a general setting in which these can be applied. 

Let pi and p 2 denote the densities for the two species, and consider a "free energy" functional T that 
depends on p\ and p 2 through 



T{pi,p2)= / F(p 1 (x),p 2 (x))dx + / / V(x -y)pi(x)p 2 (y)dxdy 

J A JAJA 



(1.1) 



The function F represents the free energy density of a system with only short range interactions. We require 
that F be separately convex in p\ and p 2 , and that the set of all (pi, p 2 ) for which F is finite, which we call 
the domain of F, is an open convex set. Finally we require that 

"~ F(p uP2 )>0 (1.2) 



dpidp 2 

for all (pi,p 2 ) in the domain of F. An example of particular interest is given by 

F(pi,p 2 ) = j [G( Pl ) + G(p 2 ) + D( Pl + p 2 ) - p lPl - p 2P2 ] (1.3) 

where (3 is the inverse temperature, p\ and p 2 are chemical potentials, G is a convex function, and D takes 
account of short range interactions between the particles which we have taken for simplicity to be species 
independent. Examples to keep in mind are: 

G(t)=t\og(t), (1.4) 

D( ., = f -t log(l - bt), if x < ft" 1 , / L5 n 
\ +00 otherwise 

for some b > 0, where b^ 1 can be thought of as the close packing density. In this case D is defined on a 
bounded interval, and the domain of F is the set < pi, p 2 < In particular the case b — 1 corresponds 
to the lattice gas with exclusion rule. 
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As for the long range interaction terms, we assume that V(x) and U (x) are non-negative radial decreasing 
functions of x e IR d with compact support and 

/ V(x)dx = a, I U(x)dx = g (1.6) 

We further define 

1= [ |a;|[-V , (x) + -I7(a;)lda; , (1.7) 
Jtr a o 

which characterizes the length scale of the interactions. We remark that it is not essential, but convenient, 
to have the same attractive interaction U for the two species. 

Notice that F is not required to be jointly convex or afhnc in p\ and p2, although this is expected on 
physical grounds, and, in any case, in general T would not be because of the interaction terms. Convexity 
of the interaction terms would require positive defmiteness of the 2x2 matrix valued function 

-U V 

V -u 

where the hat denotes the Fourier transform.. We do not make such an assumption in this paper. Instead we 
make use in our analysis of the fact, proven later, that the free energy functional that we consider satisfies 
certain rearrangement inequalities. These give a certain monotonicity to the spatial density profile of the two 
species. We study the consequences of these for the phase diagram which originates from the competition 
between the repulsive interaction terms, which would prefer to have the two species segregated as completely 
as possible, and the entropy term, which would prefer to have the densities of the two species be uniform 
over A. 

Our aim here is to study the minimizers of this functional either under constraints on the integrals of p\ 
and pi-, or without such constraints. Specifically, let ri\ and n 2 be two given positive numbers. We define 
the sets V(ni,n 2 ) and V of constrained and unconstrained pairs of densities: T>(n\,ri2) and V by 



X>(ni,n 2 ) = |(pi,P2) G IR+ x IR+ : J pidx = m and J p 2 dx = n 2 j 



(1.8) 



V= |J V( ni ,n 2 ). (1.9) 

ni,ri2>0 

The canonical (Hclmholtz) equilibrium free energy per unit volume or the grand canonical (Gibbs) equi- 
librium free energy per unit volume will be obtained by minimizing T(pi, p 2 )/\A\ under the constraint (1.8), 
or the unconstrained as in (1.9). In the former case the pi are irrelevant and we shall set them equal to zero. 

The two cases will give different equilibrium density profiles Pi(x), defined as the minimizer of J 7 in (1.1), 
when there is a coexistence of phases, i.e. when the minimizer for specified chemical potentials is not unique. 
This may happen when J-(pi,p2) is not strictly convex and the system undergoes a phase transition. 

In the case where the diameter of the torus is large compared to £, the interactions become approximately 
local. In this case, for densities p\ and p 2 that are effectively constant on the length scale I over most of A, 
3 r (pi,p 2 ) is well approximated by J r o(pi,p2) the MFT free energy functional where 

J 7 o(Pi,P2) = J^F(p 1 (x),p 2 {x))dx + a J^pi{x)p 2 {x)dx - 7} J^Pi( x ) + pl( x )] dx ■ l 1 - 10 ) 

When pi 7 p 2 are restricted to be constants, ni,n 2 respectively, then J r o(m,n2) is essentially the function 
studied by van der Waals [21] and by Korteweg [10]. (Korteweg however considered the case a < so 
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the rearrangement inequalities used here would not apply in his case.) There are interesting questions 
regarding the minimization of J-q in case pi, p 2 are not restricted to be constant. For example, knowledge 
of the minimizing profiles enables one to compute the equilibrium internal (interaction) energy at given 
temperature under the constraints in (1.8). An application of this is given in Section 6. 
We will show in section 2 that 



inf T{p\,p 2 ) and inf J r o(pi,P2) (I- 11 ) 

(pi,P2)eV (p 1 ,p 2 )G'D 

are both attained, and the same statement holds for the constrained problem. Clearly 

/o==4|, inf M P i,p 2 ) (1-12) 

|A| (pi,P2)€X> 

is independent of A. We shall see in section 3 that, using of the rearrangement inequalities, we also have 

iAi im TIT, ' m L -^(Pi'flO = fo > ( L13 ) 

A|^oo |A| {p 1 ,p 2 )GT> 

and moreover there is a close relation between the minimizers in both infimums in (1.11). 

The equality in (1.13) may look natural, but an example of Gates and Penrose [8] shows that the analogous 
statement need not hold for simple one component systems. Specifically, their one component system free 
energy functional is given by 

G(P) = ~a I P{x)^ogp{x)dx + I [ p(x)V(x - y)p(y)dxdy (1.14) 
P J a Jaja 



and 



Qo{p) = \ I p{x)logp{x)dx + a [ [ p{x) 2 dx (1.15) 

P J A JAJA 



where a and V are related as before. Since a > 0, Go is strictly convex, and hence is minimized at the 
constant density p(x) = n. However, if V(k ) < for some k , then for sufficiently small e > 0, there is a k 
close k and a S > so that with p e (x) = n + esm(kx), 

-^G(Pe) < ^o(n) - S (1.16) 

for all A sufficiently large. In our model, the rearrangement inequalities prevent oscillations from lowering 
the free energy. 

In more detail, observe that 

Hpi,pz) = Fo(pi,p2) - J J A &( x >y) v ij( x - y)^( x ^y) dxd y ( L17 ) 

where 

Vn = V 22 = -U and V12 = V 2 i = V , (1.18) 

and 

&(x,y) = Pi( x ) - Pi(y) ■ (i-i9) 
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The rearrangement inequalities will ensure that final integral in the right of (1.17) is a surface term. This 
will be proved in Section 3 below. 

The rearrangement inequalities have other interesting consequences for the phase diagram of our system. 
To be concrete, we will investigate these in the case of the van der Waals free energy functional given by 
(1.3) and (1.5). Because of (1.13), the rearrangement inequalities allow us to draw conclusions about the 
phase diagram of J-'q even though rearrangement of p\ and pi may not affect the value of J-"o(pi, pi). 

We will see from the structure of the minimizers that the system, for appropriate choices of F, V and 
U, exhibits two phase transitions: A liquid-vapor transition and a segregation transition. In the first of 
these, below a precisely determined critical temperature, the system separates into two regions, one a liquid 
region in which the total bulk density is uniformly very close to some number pi, and the other a vapor 
region in which the total density is very close to some number p v . This is similar to a van der Waals type 
transition in a one component system. The other transition is a segregation transition. The liquid or vapor 
homogeneous phases can split into two sub regions in which the densities of the two species are very close 
to certain preferred concentrations p\ and pj for the liquid, and p+ and p~ for the vapor. The values of 
these densities depend on the temperature and the parameters of the system, and a fairly detailed analysis 
is required to determine them in any given case. Nonetheless, there are general conclusions we can draw as 
a consequence of the rearrangement inequalities; for example, we will see in broad generality that for our 
system one always has 

Pj < Pi <pt< P\ (1-20) 

and never 

Pi <pi< Pi < pt ■ (i.2i) 

The paper is organized as follows: In section 2, we state precisely the rearrangement inequalities referred 
to above, and apply them to prove the existence of minimizers. We also prove that minimizers necessarily 
have certain regularity properties. In section 3 we prove a general theorem showing that as a consequence 
of the rearrangement inequalities, minimizing T and minimizing T$ are essential by the same problem for 
L much larger than £. The general theorem requires certain assumptions on F. These assumptions are 
verified for the physical model that we study in the next two sections. Section 4 is devoted to the application 
of convex analysis methods to deduce general features of the minimizers for our system. We will see that 
there are generically one, two, three or four coexisting equilibrium states in this model. The general convex 
analysis does however permit more bizarre coexistence of phases. In section 5 we impose further restrictions 
(physically natural) and obtain a detailed phase portrait, ruling out these bizarre non-generic behaviors. 
Section 6 gives an application of the results obtained here to the study of the equilibrium states of an energy 
conserving kinetic model. In particular, the results obtained here allow us to determine the equilibrium 
partition of the total energy into kinetic and interaction terms. Thus, given the energy, we are able to 
determine the equilibrium temperature in such a system. 

2. Existence and Regularity of Minimizers. 

In this section we show that minimizing densities exist, and begin the determination of their nature. 
Since particles of the same species attract and of different species repel, one can expect that the minimizing 
densities should be "well organized" in some sense, and we use a rearrangement inequality to prove this. 

Introduce coordinates x%, i = 1, . . . , d with — L < Xi < L for each i = 1, . . . , d where L > 0, and we make 
the usual identifications to parametrize our torus A. 

We say that a function g on A is symmetric monotone decreasing about the origin 0, . . . , if it is a 
symmetric decreasing function of each coordinate Xi for — L < Xi < L. That is, for each i 

g(xux 2 , ...,Xi,...x d ) = g(xi,x 2 , • • • , -x u ...xi) (2.1) 



5 



and whenever L > yi > Xi > 0, then 



g{x 1 ,x 2 , ...,yi,...x d )< g(x 1 ,x 2 , . . . , x u . . . x d ) . 



(2.2) 



We define symmetric monotone increasing in the analogous way. Finally we say that a function is symmetric 
monotone if it is either symmetric monotone decreasing or increasing. 

Given a non-negative integrable function g on A, we say that a non-negative intcgrable function g* is a 
symmetric monotone decreasing rearrangement of g if it is a symmetric monotone decreasing function such 
that for all t > {x\g(x) > t} has the same measure as {x\g*(x) > t}. Analogously we define a symmetric 
monotone increasing rearrangement h± of a non-negative integrable function h. 

We will see below that for an appropriate choice of p\ and p 2 *, J 7 (p*,p2±) < T(p\,p 2 ), with equality 
only in case p\ and p 2 are already symmetric monotone up to a common translation. The relevance of 
these definitions lies in the following lemma, which says that the "total variation", ||V#||i, of a symmetric 
monotone function g is a "surface term". (Here and in what follows, || • || p denotes the L P (A) norm.) 

Lemma 2.1 Suppose that g is a bounded and non-negative integrable function on A that is symmetric 
monotone. Then g has an integrable distributional gradient, and the total variation of g, that is ||V<?||i, 
satisfies 



Note that the right hand side is proportional to |A|( d which is what makes it a surface term. 

Proof: We suppose first that g is smooth and symmetric monotone decreasing, and derive the bound (2.3) 
under this hypothesis. Any such function g has the property that as one goes around any of the circles where 
only one coordinate varies, say Xi, one proceeds monotonically from the minimum of g on this circle, to the 
maximum, and then back down to the minimum. Since g > 0, the difference between the maximum and the 
minimum is no more than ||g||oo- Hence, since g is smooth, 



Now integrating over the remaining coordinates one obtains (2.3). 

For the general case, one may approximate g by applying the heat kernel on the torus to it. This 
approximation preserves both the positivity and the property of being symmetric monotone. Also, the heat 
evolution is continuous in both the total variation norm and the L°° norm, so (2.3) is preserved as the 
approximation is relaxed. ■ 

This lemma will play a crucial role in determining the structure of our minimizers, since, as we will see 
below, our minimizers arc symmetric monotone, up to translation on the torus, and thus, they satisfy the 
bound (2.3) as well as a-priori L°° bounds. 

Lemma 2.2 Let F be a function on LR + x LR + that satisfies (1.2) everywhere on its domain. Let U be 
a non-negative strictly monotone decreasing radial function on LR d . Suppose also that the diameter of its 
support is less than the smallest period of A. Then for all non-negative functions g and h on A, there is a 
symmetric monotone decreasing function g* and a symmetric monotone increasing h+ so that 




(2.3) 




(2.4) 




(2.5) 



J J A J J A 
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and 



[ F(g(x),h(x))dx> [ F(g*{x),K(x))dx . (2.6) 

J A J A 



Moreover, there is equality in (2.5) only in case g = g* and h = h± up to a common translation. 

These estimates (2.5) could be deduced from related rearrangement inequalities of Baernstein and Taylor 
[3], and Luttinger [14]. However, a direct proof is provided in Appendix A. The second inequality is reminis- 
cent of an inequality of Almgrcn and Lieb [2] with the opposite sense for symmetric decreasing rearrangement 
on M n in which both functions are "piled up" around the same point. The proofs though are different. 

We now turn to the existence of minimizers for the variational problem of determining 

mi{F(jH,P2) ■ (pi, P 2) eV(n u n 2 )} (2.7) 

where T>(m, n 2 ) is given by (1.8). 

The lemmas collected above enable this to be proved in a wide variety of circumstances. Theorem 2.3 
below applies to one of these that is of particular physical interest. Mathematically, it is a very direct 
consequence of the lemmas since when D is given by (1.5), there is a trivial a-priori bound on ||pi||oo and 
Halloo - namely 6 _1 . 

Without such a term to enforce a-priori L°° bounds, and without an entropy that prevents a vacuum, 
the proof is more involved, but still works. Mathematically, this context is much more interesting, if less 
physical. To keep the focus on the main physical case, we prove here the theorem only in the simple case 
required to treat the van der Waals gas and put in Appendix B the statement and the proof for the other 
case. 

The role of the rearrangement inequalities is this: Since the interaction is not positive definite, lowering 
the interaction energy may well favor oscillations in the densities p\ and p 2 , and there is no mechanism to 
prevent a minimizing sequence from oscillating more and more rapidly so that some part of the mass vanishes 
in a weak limit, with the consequence that the limit no longer belongs to T>(ni,n 2 ). This cannot happen 
here because of Lemma 2.2. 

Theorem 2.3 Let F be given by (1.3), where in particular D has the property that D(t) = oo for t > b^ 1 
for some b > 0. Then, for all n\ and n 2 and all (3, 

T{m,n 2 ,p)= inf T{ Pu p 2 ) (2.8) 

is achieved at least at one minimizing pair ( P i,p 2 ). Any such pair satisfies the Euler-Lagrange equations 
G'(p 1 ) + D'(p 1 +p 2 )+f3V*p 2 -j3U*p 1 =C 1 and G'(p 2 ) + D'{ Pl + p 2 ) + (3V * Pl - (3U * p 2 = C\ . (2.9) 
Moreover, 

||Vpi||i JVpalli < lAI^" 1 )/^- 1 , (2.10) 

and if (pi,p 2 ) is any minimizing pair, then pi = p\ and p 2 = p 2i < up to a common translation. 
Likewise, for any pi, p 2 , let 

T{pi,p 2l l3) = inf \F{pi,p 2 ) - pi l pi(x)Ax-p 2 \ p 2 {x)Ax\ . (2.11) 
{pi,pi)ev y J A Ja ) 

Then the minimizers do exist and satisfy the above conditions with C\ and C 2 in (2.9) replaced by p\, p 2 . 
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Proof: We first observe that we may assume pi,p 2 <b 1 . These a-priori bound makes this case easier. 
The functional 

Pi,P2^HPuP2) (2-12) 

is jointly continuous in the L 1 topology by the dominated convergence theorem. 

Now let (p^ ,Pz) be a minimizing sequence satisfying the constraint (1.8)i, and such that {p^\p^) is a 
symmetric monotone pair. Such a sequence exists by Lemma 2.2. By the Helly selection principle and the 
a-priori pointwise bound, this sequence is strongly compact in L 1 x L x . Passing to a convergent subsequence, 
we have our minimizing pair. 

The Euler-Lagrange equations now easily follow, and then (2.10) follows from (2.3) and the a-priori sup 
norm bound provided by D. ■ 

3. Large Volume. 

In this section we focus on the case in which L, the period of the torus A, is large compared to £, the 
length scale of the interaction defined in (1.7). 

Lemma 3.1 Let (pi, p 2 ) be any minimizer for the free energy functional. Then, using * for convolution, 

2 2 

W \V* Pj - a Pl \dx + W \U* Pj - <y P3 \dx < C\A\W d (3.1) 
j=1 Ja j=1 Ja 

for some constant C depending only on [3 and ni,n 2 . ,4s a consequence, 



2 2 

/ PiV * Pjdx - a PiPjdx + V] / pill * pjdx - u \ p % Pjdx 

• =1 J A J A - =1 J A J A 



where again, C depends only on f3 and n 1; n 2 . 
Proof: 



< C\A\ (d ' 1 ' )/d (3.2) 



Pj *V(x)- ap 3 (x) = f V(y) [ Pj (x - y) - Pj (x)] dy = - f f V(y)V P] (x - ty) ■ ydydt . (3.3) 

J A JO J A 



Hence 



J \ Pj *V(x)-a Pj (x)\dx< ^V(y)\y\dy^j ||Vpi||i , (3.4) 

and from Lemma 2.1, 

J^\p J *V(x)- a p j (x)\dx<C^J v(y)\y\dy^ lA^-^ . (3.5) 
The same argument applies to the terms involving U. ■ 



Now fix any e > 0. A simple Chebyshev argument based on the integral bounds of Lemma 3.1 shows that, 
by choosing A sufficiently large, off of a set of "surface term size" , V * pi — api ± e and U * pi — api ± e. 
More precisely, with \A\ denoting the Lebesgue measure of a set A, 

\{x | \V * pi(x) - api(x)\ > e }| < Ce~ 1 l\A\ {d ~ 1)/d , i = l,2, (3.6) 
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and likewise for U . 

Let G e be the "good" set 



f]{x\\V*p j (x)-ap j (x)\<e}\f]lf]{x\\U*p j (x)-ap j (x)\<e}\ . (3.7) 
K j=i / V=i 

Any minimizing pair of densities (pi,p 2 ) for T with m and n 2 fixed must satisfy the Euler Lagrange 
equations 

dF 

■K— {Pl,Pl) + V * p 2 - U * pi = d 

OF (3 ' 8 » 

-Z—{P1,P2) + V * pi - U * p 2 = C 2 
OP2 

for some C\ and C 2 . 

Likewise, any minimizing pair of densities (pi,p 2 ) for Tq with m and n 2 fixed must satisfy the Euler 
Lagrange equations 

dF 



dpi 
<9.F 

C 2 = (pi , P2) + api - a pi . 

OP2 

Given a minimizing pair (pi,p 2 ) for J 7 , define functions Ci(x) and C 2 {x) through 

Ci(or) := — {pi,p 2 ) +«P2 - o-pi, 
Cpi 

^ ^ 9F / 

C-2(ar) := — - (pi,p 2 ) + cepi - ap 2 . 

OP2 



(3.9) 



(3.10) 



It follows from the definitions that on G e 

\Ci{x) - d| < e and |C 2 (a;) - C 2 | < e . (3.11) 

Therefore, on the set G £ , (pi(x), p 2 (x)) is a solution of (3.10) with C\{x) and C 2 (x) very close to C\ and 
C 2 . We wish to conclude that on G e , the values of (pi(x), p 2 (a;)) are essentially those of a minimizer of i*b. 

Toward this end, we make a "stability" assumption on the Euler Lagrange equations describing minimizers 
of T under the constraint for fixed n\ and n 2 . This condition is easy to verify for many particular choices 
of F, as we shall see. 

Definition: (Amenable Free Energy Function) We say that the free energy function F(pi, p 2 ) + apip 2 — 
(p\ + p 2 ) is amenable when for any given constants C\ and C 2 , there are at most a finite number N of pairs 
of numbers 

i = 1, . . . ,N (3.12) 

so that for any given e > 0, there is a S > depending only on e, C\ and C 2 such that the following is true: 
Whenever any pair of numbers (pi,p 2 ) satisfies 

dF dF 

-^—{pi,P2) + ap 2 -ap 1 =Ci and -— (pi, p 2 ) + a.p\ - op 2 = C 2 . (3.13) 
Cpi Cp 2 

9 



for some C\ and C 2 with \C\ — C\ \ + \C 2 — C 2 1 < S, it follows that 

\Pi ~ Pi } \ + P { 2 I < e ; for some l<i<N . (3.14) 

Checking this in practice amounts to checking that there is continuous dependence of solutions to (3.9) 
from (Ci, C2). Just to concrete, consider the simple case in which F(pi,p 2 ) = (pi log pi + p 2 logp 2 )/(3. Then 
(3.9) are 

logpi + /3ap 2 = C\ and logp 2 + /?api = C 2 (3.15) 

Then with 

hc lt cM = e Cl exp(-/foe C2 (exp(-/3ap))) , 
pi and p 2 satisfy the fixed point equations 

P = h Cl ,c 2 (p) and P = h c 2 ,c\(p) 

respectively. For each C\ and C 2 there is a number i? depending only on C\ and C2 so that hc 1 ,c 2 (p) is 
convex for p < R, and concave for p > R. One sees that there are always one, two or three solutions to (3.15) 
in this case. As C\ and C 2 vary in small intervals, there are at most three small pairs intervals required to 
hold all of the solutions. This is the situation described in general by the definition. 

We will give further physical examples in the next section, and here we focus on the general consequences 
of amenability. 

The main consequence is that on the "good set" G e , any minimizers p\ and p 2 of T are essentially discrete, 
taking their values in the union of a finite number of short intervals. We now identify these intervals, and 
relate the minimization of T and Tq. 

Theorem 3.2 Suppose that the free energy function is amenable. For any fixed (n\,n 2 ), let (pi,p 2 ) be a 
minimizer for T in T>(ni,n 2 ). Then there is a finite set of pairs of numbers 

(p^pW) i = l,...,N (3.16) 

such that for all e > 0, and all L sufficiently large, there is a set G C A such that 

|xf<qA|-v- 

such that for all x in G e , 

\pi(x)-p^\ + \p2(x)-p^\<e (3.17) 

for some 1 < i < N. Moreover, there is a pair (pi,p 2 ) satisfying \pi — p^ | + \p 2 — p^ 2 \ < e that satisfies the 
Euler-Lagrange equation (3.9) for a minimizer of in T>(ni,n 2 ) for some n\ and n 2 with \n\ — m\ + \n 2 — 
n 2 \ < e. 

Proof: For some k to be determined later, let G = G K . Then, as a consequence of Lemma 3.1, pi(x) and 
p 2 (x) satisfy (3.10) for Ci(x) and C 2 (x) nearly constant on G. Let C\ and C 2 be the respective average 
values of C\(x) and C 2 (x) on A. It follows that on G, p\ and p 2 have values lying finite number of pairs 
of intervals. As k is decreased, so is the width of these intervals. Each pair of values in these intervals is a 
solution of the Euler-Lagrange equations (3.9) for some values of C\ and C 2 close to C\ and C 2 respectively. 
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By decreasing n, we can ensure that \C\ — C\ \ + IC2 — C2 < 5, where <5 is related to e as in the definition of 
amenability. ■ 

We remark that in specific cases, it is possible to carry the analysis further, and to show that the pairs 
of values (pi,p2) correspond to values of minimizers for JT , and not simply solutions of the Euler-Lagrange 
equation. At this level of generality, that is not possible. 

In this section we have only considered the constrained problem because in the unconstrained case the 
situation is much simpler. In the unconstrained case surface tension would discourage partitioning A among 
different minimizing phases. 



4. Local Interactions. 

In this section we consider the case in which V(x) — aS(x) and U(x) = <rS(x). This will enable us to 
obtain a very complete picture of the minimizers when a is sufficiently small and a rather clear understanding 
for arbitrary a, as we show in the next section. Note that this corresponds exactly to replacing T by J-q. 

We begin our analysis by looking for the unconstrained spatially homogeneous minimizers of Tq with 
additional chemical potentials p\ and p 2 (grand- canonical ensemble) which will play the role of Lagrange 
multipliers later on. Therefore, fixing pi and p 2l we look at the minimizers of grand canonical free energy 
density on 1R + x M + given by 

/mi,M2(Pi>P2) = F {PuP2) + a Pl p 2 - ^(pl + pi) - pipi - P2P2 (4.1) 

Since this model exhibits both condensation-evaporation transition and segregation transition the natural 
variables are not pi and P2 but p and <j) where 

p(x) = Pl (x) + p 2 (x) and <f{x) = ~ (4 2) 

p(x) 

Then the free energy J-q(p\, p 2 ) can be written as 

Fo{pi,Pi) = I 9n,h{p(x), <f>(x))dx , (4.3) 

J A 

where 

g^ h (p, 0) = F (|(1 + 0), |(1 - 0)) + |p 2 (l - 2 ) - a -p\\ + <f) -pp- hp<p (4.4) 

where 

M=^, h=^. (4.5) 

It turns out that, under our assumptions, minimizers of this functional have a very special structure: 
Each of pi and p 2 can take on at most four values. The region A is decomposed into at most four subregions 
in which both pi and p2 arc constant. The four possible values of the densities result from a possible 
condensation-evaporation transition, together with a possible segregation transition. 

Let g^ t h be the convex minorant of g^^h- That is, 

9n,h(P> 4>) = sup{^(p, 4>) I t{p, <j)) < g^,h(p, 4>)} (4.6) 
eec 

where C is the set of all the linear functions £(p, <f>) = ap + b(f> + c, a, b, c e M. Let B be the set of p and <f> 
for which 

g^,h{p,4>) > g^hiP^) ■ (4-7) 
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B is the set where there are "flat spots" in the graph of g^h, and as is clear and well known, these are 
irrelevant to the minimization problem so that 



inf Jr °^ 1 ' /3 ^ 



P1,P2 |A| 



in fu7 / 9^.h(p( x )^( x )) dx ■ (4-8) 
P><t> |A| J A 



Next, define 



and 



0e[-i,i] 



i>»,h(p)= i r nf ,9p,,h{P><f>)- ( 4 - 10 ) 
<pe[-i,ij 



Then is the convex minorant of Vv.fr- Indeed, the epigraph of g^y, i.e., the set 

{(p,<j>,z) : z>g^ h { Pl 4>) } (4.11) 

is convex, and the epigraph of Vy h is simply the projection of the epigraph of g^^ onto the p, z plane. 
It is clear that 

J g„, h (p( x )^( x )) dx > J ^,h{p( x )) dx (4.12) 

with equality if and only if for almost every x, <p(x) minimizes gn,h{p( x )> 4>) considered as a function of (f>. 
Now let p avo be given by 

Pavo = |X| j A p ^ dx ■ ( 4 - 13 ) 

Then by Jensen's inequality, 

j^l J^v,h(P( X )) dx > ^,/i(Pave) ■ (4.14) 

It is well known that if ipp,n,h is strictly convex, then there is equality if and only if p(x) = p ave almost 
everywhere. However, even if ip^h is not strictly convex, the proof of Jensen's inequality has important 
consequences for our problem. 

To explain these, we first recall that if tp is any convex function on the positive axis, then an affine function 
£(p) = ap + b is a supporting line for ip in case ip(p) > £(p) for all p, and if ip(po) = £{po) for some po- If £ is 
any supporting line for tp, then the set 

{p : £(p) = V(P) } (4-15) 

is a closed interval. Such an interval is called a support interval of ip. The following is simply Jensen's 
inequality, with the only novel feature being that the statement about the cases of equality that applies 
outside the strictly convex case. 

Lemma 4.1 Let (fi, <S, z/) be a probability measure space, ip a convex function on the positive axis, and p a 
non-negative measurable function. Then 



tp{p)dv >^yj P d ") . (4-16) 

and there is equality if and only if, up to a set of measure zero, p takes it values in a single support interval 
oftp- 
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Proof: This follows from a close examination of the standard proof of Jensen's inequality, which turns on 
the fact that if / and g are two measurable functions, 



J max{/, g}Av > max j J fdv , J gdv j , (4.17) 

with equality if and only if either / = max{/, g} or g = max{/, g} almost everywhere. This applied to 

</>(/>)= sup (4.18) 
e 

where the supremum is taken over all supporting lines of ip yields the inequality. By the above, there is 
equality if and only if for any two supporting lines £\ and £ 2 , cither £i(p(x)) = max{£i(p(x)) , £2(p(x))} or 
^2{p(x)) = max{£i(p(x)),£2{p(x))} for almost every x. This can only happen if, almost everywhere, p takes 
on all of its values in a single supporting interval. ■ 

This has the immediate consequence that for any minimizcr of the local interaction problem, the values 
of p{x) all lie in a single support interval of Vv,/i- Now in most cases that we will consider, 

MP)>^AP) (4-19) 

for all p in any support interval, except at the endpoints. In particular, this is the case if il>n,h is almost 
everywhere either strictly convex or strictly concave, and there are no lines that are tangent at three or more 
points. Then a minimizer for our problem cannot have p(x) in the interior of the support interval because of 
(4.19). In exceptional cases, there may be physical values in the interior of a support interval. This occurs 
when there is a point of triple tangency, or higher. 

Putting aside this exceptional case for the moment, there are only two possibilities for any minimizcr in 
our problem: Either 

• The support interval of ip^h that contains p avc consists of p avc alone. In this case, p(x) = p avc almost 
everywhere. 

• The support interval ofip^,h that contains p avo consists of a closed interval [p v ,pi] with p v < pi. In this 
case, A = A„ U A{ with p(x) = p v almost everywhere in A v and p{x) = pi almost everywhere in A^. 

We say that there is a condensation-evaporation transition in the second case; the region A„ contains the 
vapor state and Ai contains the liquid state. The volumes of these two regions are given by 

\At\pt + \Av\pv = |A|p avc = m + n 2 • (4.20) 

Note that the volume fractions of the liquid and vapor states arc determined by this relation alone, before 
we begin considering any possible segregation in cither the vapor or the liquid. 

We note that there is a condensation-evaporation transition exactly where Vv>fr nas an interval of concavity. 
Any such interval is contained in an interval [p v ,pt\ with 

ipit,h(Pv) = ^^,h(Pv) and 4>n,h(pi) = $n,h{Pt) ( 4 - 21 ) 

and 

This means that 



pi 

< h dp - (4.22) 
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which is Maxwell's equal area rule for determining the values of p v and pg. 
Suppose for all t e (0, 1) 

^M 1 - *)Pv + tPt) > (1 - t)^M(pv) + t^,h(Pe)- (4.23) 

Suppose also that the infimum in (4.9) is attained exactly at one or two values of <f> when p = p v and p = pg. 
Then, depending whether p v equals or not pi and whether there are one or two minimizers for <j> in (4.9) there 
will be one, two, three or four states. To determine if these hypotheses hold we must examine a particular 
free energy. 



5. van der Waals Gas. 



In this section, we focus on the van der Waals gas for which 

fp,^ 2 {pi,P2) = ^[G(p!) + G(p 2 ) + D(p)] +ap lP2 - |(p? + pi) - PiPi - P2P2 (5.1) 

and hence 

9fi,*,h{p, <t>) = -p\P (f (! + ^))+G (f (1 - 0)) + D(p)] + -^p\\ <j?) °-p\\ + <?) -pp- hp<i> (5.2) 

The function G, defined in M+, is assumed to be smooth and moreover 

(1) G and G" strictly convex; 

(2) G'(x)->-oo, G"{x) +oo for x 0. 

The function D is defined in some interval (0,6 -1 ) contained in (and possibly coinciding with) M + , and 
is smooth and strictly convex. For example, with G(t) = t\og t and D given by (1.5), these conditions are 
satisfied. This is the usual two components van der Waals gas. 

The existence of minimizers follows from the considerations of Section 2, which apply also in the case of 
local interactions. Note however that we have no symmetric monotonicity of the minimizers and, in general, 
no reason to have regularity properties. Moreover, the same argument that we used in the proof of Theorem 
2.3 shows that any minimizer satisfies the Euler-Lagrange equations 

G'(p 1 )+D'{p 1 +p 2 )+[3ap 2 -f3<jp 1 =C 1 
G'{p 2 ) + D'( Pl + p 2 ) + /fapi - pa/* = C 2 . 

We now turn to the study of the structure of these minimizers. The following theorem gives a criterion 
for segregation: 

Theorem 5.1. With g given by (5.2), the infimum in (4-9) is attained in exactly one point <j>(p; (3, p, h) if 
h ^ 0, independently of p, or if h = and [3(a + a) < G"(p/2). If h = and (3(a + a) > G"(p/2), then the 
infimum is attained in two points <p± — ±(f)(p; (3, p, 0) 

Proof: When there is no ambiguity, we omit subscript [3, p, h from gp^,h- For fixed values of (3, p and h we 
look for the minimizers of g. To do this, we first fix p and look at the solutions of the equation 

|?GM) = (5.4) 
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for fixed p. Let <p — <fi(p\ P, p,h) be one of the solutions to (5.4). Explicitly, for fixed p > 0, (f> solves the 
equation 

H p (4>)=[3{a + <T)p4) + 2h, (5.5) 

with 

W): =G'(|(l + 0))-G'(|(l-^)) (5.6) 
The following properties of the function H p will be relevant below. Assume p > 0: 

H p {-4>) = -H p (4>), H p (0)=0- (5.7) 

H' p {<f>) >0,V</>e (-1,1); (5.8) 

H p ((t)) -> ±oo, i?p(0) -> +oo for <j) -> ±1; (5.9) 

0i/;(0)>O for0#O, H;(0) = 0; (5.10) 
(Notice that if G(£) = ilogi then H p {4>) — 2 tanh -1 (</>)). In particular, (5.10) implies that 

inf H' p (4>)=H' p {0)=pG"&) and H p (d>)>pG"£) if <P^0. (5.11) 

<f>e[-i,i] F 2 2 

The above relations are immediate consequences of the assumptions on G: in fact (5.7) just follows from 
(5.6); (5.8) follows from 

H'M) = § G "(f (! + *)) + G "(f (! - *)) ( 5 - 12 ) 
and the fact that G is strictly convex. (5.9) follow from condition (2) on G and finally, (5.10) follows from 

H';W = ^G'"(f (1 + 0)) - G'"(|(l - 0)) (5.13) 

and the strict monotonicity of G'". 
We have 

0(p, 0) = | - (a + CT )/3p) > ^ (g"(|) - (a + a)/?) . (5.14) 

Therefore, if 

f3{a + a)<G"{ P -), (5.15) 

the r.h.s. of (5.14) is positive, the function g(p, <p), for p fixed is a strictly convex function of <j) and hence it 
has a unique minimizer solving <p(p; j3, p, h) solving (5.5), because of condition (2) on G permits to exclude 
that it is a corner solution. On the other hand, if 

/3(a + CT )>G"(|), (5.16) 

by (5.10), there is an interval (—(f) s ,<f> s ) where the function g(p,<j>), for p fixed is concave while in the 
complement it is convex. The value (f) s is determined by the condition 

H' p {<j> s ) = {a + <j)(3p. (5.17) 
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Therefore, it is possible to get more than one stationary point, satisfying (5.5). 

We assume that (5.16) is verified and look at the stationary points. We distinguish two cases: h = and 
h ± 0. 

We assume first h — 0. Then, clearly, <j> = solves (5.5). The condition (5.16) ensures that g{p,4>) is 
concave in <j> — and hence <fi = is not a minimizer. Moreover, by the symmetry, if (p(p; [3, p, 0) solves 
(5.5) the same is true for — <f>(p; (3, p, 0). (5.8) and (5.10) show that there is exactly one positive value 
(j>(p; (3, p, 0) € (<j) s , 1) solving (5.5). Consequently, we have the two stationary points 

<j>± = ±$(jr,0,H,O). (5.18) 

Now we take h ^ 0. One immediately realizes that there is h(p) such that, if h £ [-h(p), h(p)] then there 
is only one solution to (5.5), for h = ±h{p) there are two solutions, one of them in (— <f> a , <f) s ) and the other one 
is the minimizer. Finally, if h e (—h(p),h(p)) there are three solutions, one in (—<t> 8 ,<f> a ) an d the others, in 
the complement, corresponding to local minimizcrs; moreover they have different signs. Let 4>i(p; (3, p, h) > 
and §i{p\ (3, p, h) < denote the two local minimizers. We show that g(p, <j>i{p; (3, p, h)) ^ g(p, (faip; (3, p, hj) 
and hence there is only one absolute minimizer. To show this, we compare with the case h — 0, where 
<j>2(p; (3, p, 0)) = — (f>i(p; (3, p, 0)) and the corresponding values of g are equal. We define 

JAh)=gp, ll „(p,<j>i(p-,P,ii,h)), i = l,2. (5.19) 

Since are stationary points for g, we have 

A. j, _ S&« ta Ur, ft * k))% + ft ft. ,, »)) (s 

= ~^f^(' ' A m, M) = -pMp; P, p, h). 

Therefore, for h > 

Ji(ft) < Ji(0) = J 2 (0) < J 2 (/i) (5.21) 

and vice-versa for /i < 0. ■ 



Next we give a criterion for condensation-evaporation transition. Such a transition occurs exactly for those 
values of a, a, (3, p, h for which the function ipp lllt h(p) fails to be strictly convex. Our first result pertaining 
to this is the following: 

Lemma 5.2 Let 4>{p) be any minimizer of g{p,4>) in (5.2) with respect to <f>, so that ip(p) = g{p,<f>{p)). 
Then ijj"(p) is strictly positive if the point (p, (p(p)) is in the regions where the Hessian of g, D 2 g, is positive 
definite. 



Proof. 

We have: 



dp 2 d<f>dp dc/) 2 d<p ^ 

= ((!,$), D 2 g(l,$)) + ^", 
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D 2 g being the Hessian matrix of g and ( • , • } the usual scalar product in IE 2 . 
When 6{p) is a minimizcr, above relation reduces to 

iP" = ((l,$),D 2 g(l,2)) (5.24) 
Therefore, ip"(p) is strictly positive if the point (p,6(p)) is in the regions where D 2 g is positive definite. ■ 



The above result can be restated in terms of the function /. Since 6(p) solves (5.5), it is plain that 
V{p) = <J T (W), D 2 f j T {l,h) where 





/ dp 


dp \ 




3 = 


dpi 

86 


8p 2 


(5.25) 




\8p~i 


8pi ) 





Therefore we study the positivity properties of D 2 f. Its expression is: 

D 2 f _ (G"( Pl )-f3a + D"(p) D"(p) +/ 3a \ 

f -{ D"(p) + (3a G"( P2 ) -(3* + D"(p) J ■ ^ 2b) 

As a consequence of Lemma 5.1, when the Hessian of / is positive definite, there is no condensation- 
evaporation transition. However, if a is sufficiently large, the diagonal terms may become negative, while 
for a sufficiently large the off diagonal terms may make det(D 2 /) < 0. 

We now consider explicitly some examples for the case G(t) = t log t. 
1) a = and D = (ideal gas). Then D 2 f is positive definite in the region V defined by 

P = {(puP2)\piP2 <a- 2 p- 2 }. (5.27) 
To check this, first suppose aj3p > 2. With our choice of G, (5.5) reduces to 

= tanh (^P>) . (5.28) 

With £ = a[3p/2 it becomes 

6 = tanh(£0) . (5.29) 

For all strictly positive values of 

sinh(^) > £,6 . (5.30) 

By (5.29) we can divide the left hand side by tanh(^) and the right hand side by 6 without affecting the 
inequality. It follows that cosh(^) > £, and hence that (1 — tanh 2 (^0)) _1 > £ 2 . By (5.29), this means 
that 

(W 2 )- 1 ^) 2 (5.31) 

whenever a[3p > 2. Of course if a[3p < 2, then, by Theorem 5.1 6 — 0, and so (5.31) holds in all cases, 
and there is equality if and only if a(3p = 2. Since 

PiP 2 = ^(l-0 2 )<^, (5.32) 
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by (5.29), the condition in (5.27) is fulfilled for all the values of the parameters but for a(3p = 2. In 
this case however, = and by (5.23) ip" = - > 0. Hence 4>{p) is strictly convex and there is no 
condensation-evaporation transition for any value of the parameters. 

2) D{p) is given by (1.5) and a — 0. Then the strict convexity of G and D ensures that the diagonal terms 
are positive and we only need to check det(D 2 /) > 0. Hence, in this case D 2 f is positive definite in 

Q = {(pi,p 2 ) | det(D 2 f) >0} (5.33) 

and its boundary, 

dQ={(p 1 ,p 2 )\dct(D 2 f) = 0} 

separates it from the region where D 2 f is not positive definite. Explicitly, the condition dct(Z) 2 /) > 
can be written 

— -a 2 (3 2 +D"{p){— + — -2a/?) > 0. (5.34) 

PlP2 Pi P2 

Note that the set Q contains the set 

V = {(pi,p 2 )\ — -a 2 (3 2 >0} (5.35) 
P1P2 

because the arithmetic mean of two positive numbers is not less than the geometric mean. Since the 
equation for <f) is unaffected by the presence of D, ip(p) is strictly convex and we get the same conclusion 
as before. 

3) If a > 0, the set V is replaced by 

V a = {(pi, pa) I (- - M(~ - M - u 2 (3 2 > 0}. (5.36) 

Pi P2 

Therefore, as a increases the set V a shrinks thus encouraging the condensation-evaporation transition. 

Korteweg [10] has also discussed a situation with segregation and four phases in equilibrium in a mixture 
of two van der Waals gases. Korteweg's paper concerns the a. < case where there is attraction between 
different species too, and to which rearrangement arguments do not apply, though segregation is still possible 
because it is controlled by a + a. 

Up to now we just gave general conditions for the absence of condensation-evaporation transition. We 
now derive a formula for the pressure which allows us to draw conclusions about the not strict convexity of 
the function ip. Introduce the specific volume v = p' 1 and define the function 

q(v) = W>(v _1 ), (5.37) 

which has the same convexity properties of ip in the corresponding points. A simple calculation shows that, 

q(v) = log(v -b)- + 7(£(»- 1 )) (5-38) 

with 

7 W - \ log ^-Pj + \ log (i±|) - H P (5.39) 

and 

= a - a + 4> 2 (a + a). (5.40) 
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Now, by (5.5), we have 



^ + 7 ' = 0. (5.41) 



Hence the pressure p as function of the specific volume is given by 

«„):=-,'(„) .-L--*&pl. (5 . 42 ) 

This expression is similar to that for the one component van der Waals gas, but for the dependence of a on 
v through (j). 

Several cases are possible, depending on a and a. Consider for example the case a. > a. Then, in 
the absence of segregation, the pressure is a strictly decreasing function of v. Since the pressure must be 
equal in the liquid and vapor phases, this precludes a condensation-evaporation transition in the absence of 
segregation. (This is quite different from the usual van der Waals case, or the case considered by Korteweg) . 
The segregation permits a condensation-evaporation transition, and lowers the pressure in the liquid phase 
from what it would be without segregation. As the temperature is lowered further, the vapor phase may 
also undergo segregation. In this situation, the volume fractions are not determined. However, in the final 
section of the paper we shall see that because the local interaction model is the large volume limit of a model 
in which the minimizers must be symmetric decreasing, we can say more about this case. Thus when a > a, 
we have three critical inverse temperatures (3\ < (3^ < /?3- The first corresponds to a segregation transition. 
At the second there is a condensation-evaporation transition at which a segregated liquid phase is produced, 
and an unsegregated vapor phase. At the third, the vapor phase as well segregates. As a consequence of the 
above formula for the pressure, we prove the following: 

Theorem 5.3. Suppose that a < a . Then there is a critical inverse temperature (3 C such that for (3 > [3 C a 
condensation-evaporation transition occurs, with or without segregation, depending on the values of (3, h and 
mu, according to Theorem 5.1. On the other hand, if a> a, then either there is a unique minimizer or the 
condensation- evaporation happens at lower temperature than segregation. 

Proof: First wc consider the case a < a. We differentiate (5.42) and get 

PP («) = -(^5 + + • ( 5 - 43 ) 

Since <j> is non decreasing, the last two terms are both positive and at least linearly increasing with (3. Hence, 
for any fixed v and sufficiently large (3 they dominate the first term. On the other hand, for any fixed [3 and 
v sufficiently close to b, the first term dominates. Therefore the pressure is not monotone. 

Let us now consider the case a > a. In the absence of segregation, the pressure is a strictly decreasing 
function of v and this precludes a condensation-evaporation transition in the absence of segregation. (This 
is quite different from the usual van der Waals case, or the case considered by Korteweg). The segregation 
permits a condensation-evaporation transition, and lowers the pressure in the liquid phase from what it 
would be without segregation. As the temperature is lowered further, the vapor phase may also undergo 
segregation. ■ 

In conclusion we may have at most four phases, characterized by the following values of the density: 

Pv = PV 7y , PT = Pi O • ( 5 ' 44 ) 
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We now return to the search for the constrained minimizers which are density profiles p\(x), p 2 (x) such 
that 

/ dxpx{x) = ni|A|, / dxp 2 (x) = n 2 \A\. (5.45) 

J A J A 

Clearly if min{ni,n 2 } < p~ or max{iii,n 2 } > p\ the only possible solutions are p\ — n\ and p 2 = n 2 . 
Otherwise, non homogeneous minimizers are possible. Because of the previous lemma it is clear that p\ 
and p 2 can only take the above values in four regions A, B, C and D whose volumes are the only relevant 
properties because we are dealing with the local interaction case. The minimizer (p\, p 2 ) is given by 

!pf if x e A, ( p~ Mx&A, 

p~ if x e D, { p+ if x e D. 

or viceversa. Let a, 6, c and d denote the ratios of above volumes with |A|. The constraints are then 

api~ + bpl + cp+ + dp~ = m, apl + bpf + cp~ + dp+ = n 2 , (5-47) 
which, together with the relation 

a + b + c + d=l (5.48) 

do not suffice to determine the fractions occupied by each phases. Note that, in the case (3 s (pi) < /3 < (3 s (Pv) 
we have p+ = p~ and the relations available are sufficient to determine the volume fractions occupied by the 
three phases. 

Now we return to the finite volume case where the rearrangement inequalities take over and ensure 
regularity of the minimizers. In fact because the minimizers must be symmetric monotone, not all of the 
phases can be in contact, and there are constraints on the ordering of the densities. Suppose for example 
that we have 

p+ > pj > P + > p- . (5.49) 

We shall show that this is impossible. 

Pick numbers a, b and c separating these density values so that 

p+ > a > p~ > b> p+ > c> p- . (5.50) 

Define 

A = { x : pi(x) > a } B = { x : a > p\ (x) > b } 
C = { x :b> px{x) > c } D = { x : c> pi(x) } . 

By the rearrangement inequalities, we see that p 2 must take its minimum in A, but we know from (5.46) 
that p 2 must take the value pj in A. So we must have instead that 

Pi > Pv > Pv > Pi ■ (5-51) 

It is not evident how to deduce this in such generality by direct consideration of the local interaction model. 
The above regions are disjoint. To further determine the structure we need to consider the surface tension 
across their boundaries. In a later paper we will rigorously investigate the surface tension. Here we proceed 
on the assumption that because of it the phase boundaries will be so arranged as to have minimum surface 
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area. Then, for a two dimensional torus, these regions would be arranged in the manner indicated in the 
diagram below: 



Note that on a torus, the solutions of the isoperimetric problem are either disks, or bands depending on 
the size of the area to be enclosed. 

Now we see that the liquid state rich in species one can only be in contact with the vapor state rich 
in species one, which in turn can only be in contact with the species one rich liquid state and the species 
two poor vapor state. In the case that all four phases are present, it is the surface tension of the domain 
boundaries that fixes the volume fractions in the finite range case. 

6. The Microcanonical Minimization Problem. 

In the previous sections we have considered a minimization problem related to a model in thermal contact 
with a reservoir. In this section we extend our results to a kinetic model preserving energy, but, for sake 
of simplicity, we confine ourselves to the ideal gas case with repulsion between different species and no 
attractive part. For example, consider the Boltzmann-Vlasov equations in [4], for the position and velocity 
distributions fi(x,v,t). The results on the minimization problem we have considered in the previous sections 
can be applied to determine the equilibrium solutions of these equations. 

The equilibrium solutions of such a model are the minimizers of the entropy with energy and mass 
constraint. Therefore, the variational problem we want to deal with is the following: find the minimizers of 




c 



B 



C 




(6.1) 



in the set 




(6.2) 



where /i and /2 are probability densities in the phase space A x IR 3 and 




(6.3) 



21 



This is equivalent to minimizing without constraints the functional 



2 

0(f) = S(f) + 0E(f) -E^./*/ a dvft 



(6.4) 



on 



■D={f=(f 1 ,f 2 )\\f i €L+(AxIR 3 )} 



(6.5) 



with (3, Hi, Lagrange multipliers to be determined by the constraints. 

As is well known, the Euler-Lagrange equations for this functional force the w-dependence of fi to be a 
Maxwellian at inverse temperature fi > and hence we arc reduced to the following variational principle: 



2 

F(P\,Pi)= / drr Vpi[logp, -JXi] + (3 / dxp 1 V*p 2 , 
Ja i=1 Ua 



on 



V = {(p 1 ,p 2 )\\p i £L+(A)}, 
with [3 > 0, p,i Lagrange multipliers to be determined by the constraints 



E( Pl ,p 2 ) = \A\ 1 J A dx ^P 1 (Pi+P2) + J dxp 1 V*p 2 
Here fli — pi + log(y / 2~7r/iF r '~) 3 . This is equivalent to minimizing 



= e, 



|A| 1 / dxp t = m 
J A 



S{PUP2) 



dxp l log pi 



on the set 



{(p 1 ,p 2 )\\M~ 1 dxpi=ni, i = 1,2; E(f) = e}. 



(6.6) 



(6.7) 



(6.8) 



(6.9) 



(6.10) 



The discussion in previous sections obviously extends to the study of the minimizcrs of the functional 
(6.6). Therefore we now only deal with the solvability of the conditions on the Lagrange multipliers. This 
relies essentially on the local interaction case, because the arguments of section 3 apply. 

Consider first the local interaction case, with unconstrained energy. For [3 sufficiently small, (a[3n < 2) 
there is only a homogeneous solution. The energy as a function of fi is given by 



E(J3) = -(i-'n + an 2 \ A|, 



(6.11) 



i.e. the usual linear behavior in [3 1 — T. For sufficiently large [3 (a[3n > 2) the solution is of the form: 



Pi(x) 



x e A, 



P2(X) 



p- x £ A - A, 
with p± = p(l ± (j>((3))/2, 4>([3) the unique positive solution to 



p_, x £ A, 
p+ x £ A - A, 



(6.12) 



tanh- 1 m = ^(f3), 



(6.13) 
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and A a suitable subset of |A|, whose volume is determined to fulfill the mass constraints. Of course in the 
present situation p = n = m + n 2 and hence the energy is given by 

E(P) = + aj^dxp lP2 = l^n + a\A\ P+P _ = + ^!(1 - ^f). (6.14) 

This is obviously continuous because </>(/3) — ► as [3 — > 2/na. Moreover by using arguments similar to those 
of section 5 one can show that 0(/3) is increasing for (3 > 2/na. Therefore E{[}) is monotone decreasing, as 
sum of two decreasing functions, and hence invertiblc. 

The function /3(e) is well defined and, with this definition, the phase space densities for the two species 
are given by 

e -/3(e)„ 2 /2 

where pi are the minimizers of the free energy at temperature /3(e). 

These densities give us the minimizers of the entropy S(f) under the constraints E(f) = e. Note that 
the determination of the function /3(e) in (6.15) require knowledge of the partition of e into its kinetic and 
interaction parts. This is provided by the results in the previous sections. 
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Appendix A. 

Here we prove the rearrangement inequalities used in the paper. 

Proof of Lemma 2.2 We shall apply a sequence of symmetrization operations, and show that unless the 
minimizer is symmetric monotone, then one of the these operations would strictly lower the integrals in (2.5) 
and (2.6). 

These symmetrization operations are "transplants" to the torus of the following symmetrization operation 
on the circle S 1 = {(£1,2 ) | x\ + x\ = 1} in the X\,X2 plane. Fix any unit vector u on S 1 , and define the 
reflection R u on S 1 given by 

R u x = x — 2(u ■ x)u (A.l) 

where x — (x,y), and the dot denotes the usual inner product. It is also convenient to define R u 9 for 
— 7r < 9 < 7r by 

#„(cos(0),sin(0)) = {cos(R u 9),sm{R u 6)) . (A.2) 
Next, define an operator, denoted R+, on measurable functions on S 1 by 

p+ a ( x ) = I ma - x {a{x), g(Rux)} if x ■ u > , (A , 

u9K ' \min{g(x),g(R u x)} if x ■ u < K 

where uq denotes (0, 1). That is, the line through the origin perpendicular to u divides the circle in two, and 
the two halves are identified by the reflection fixing this line. The symmetrization operation swaps values 
at reflected pairs of points, if necessary, so that the large value is always on the side containing the "north 
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pole", (0, 1). We also define an operator R~ in the same manner, except that we put the small values on 
the side with the "north pole" . In other words, we interchange the minimum and maximum in (A. 3) 
We now state the lemma for which these definitions were made. 

Lemma A.l Let F be a symmetric function on M + x IR + that satisfies (1.2) everywhere on its domain. 
Let K be any strictly decreasing non-negative function on IR + . Then for any two bounded non-negative 
measurable functions g and h on S 1 identified with [—n, n), 

r r g(8)K(\e - mmed^ > f f R+ g (e)K(\e - <j>\)Rzh{4>))d9d<t> (a.4) 

J— TT J — TT J —TT J —TT 



and 



r F(g(9), h(8)d9 > r F(R+g(9), R-fi(6))d9 . (A.5) 

J— TT J —TT 



Moreover, there is equality in inequality (A. 4) if and only if for some fixed 6q, R^g{0 — 9 ) = g{0) and 
R~h(6 - 6>o ) = h(6) for almost all 9. 

Lemma A.l will be proved after using it to prove Lemma 2.2. The argument is adapted from [3], which 
however does not consider cases of equality. 

Consider first (2.5). Fix any index i, and fix values of Xj and yj for j ^ i. Determine a sequence of unit 
vectors {uj} inductively as follows. Suppose that the Ui have been determined for i < j. For i < j, define 
gi and hi inductively by 

9i = Rti9i-i, 5o =.9 and h i = R~.h i - 1 , h = h . (A.6) 

Now choose Uj so that 

J —TT J —TT J —TT J —TT 

~\* U tF{{[ f R u9(6)K(\9-4>\)RZh(4>))dedct>- j V ^ g 3 ^(9)K{\Q - 4>\)hj-iim0^ 

That is, we choose the direction uj to get an effect that is at least half of the largest possible effect at that 
stage. 

The operations R^. preserve the modulus of continuity, so if g and h are continuous, the sequences {gj} 
and {hj} are strongly compact. Passing to a subsequence along which limits exist, define 

g* = lim gj and h± — lim hj . (A. 7) 

j— >oo j^oo 

It follows from the choice of the sequence that for every u, 

r r R+g*(e)K(\9 - ^R-.h+mdow - r r g^emio - m*(modct> . (a.8) 

J —TT J —TT J —TT J —TT 

By Lemma 5.3, this is only possible if g* is symmetric monotone about the "north pole" , and h± is symmetric 
monotone about the south pole, — u . Thus, since the {gj} are equimeasurable, as are the {hj}, g* and h+ 
are the symmetric monotone rearrangements of g and h respectively. Since 

f T R+ gj -i(9)K(\6 - MRzfr-^MMM - r g^mx^ - ^hj-^dedep < (A.9) 

J —TT J —TT J —TT J —TT 
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for each j, it follows that 



r r g*(9)K(\e - ^)/i*(0))d0<ty < r r <?(#)# - . (a.io) 

J —TV J — IX J— TV J — 7T 

The requirement of continuity is easily removed by a density argument, exactly as in [3]. This proves the 
inequality in general. Now for equality to hold, it must be the case that 

^ r K9(0)K(\9 - c/>\)R-.h(cf>))d6dct> = r r g{6)K{\6 - <j>\)h(4>))Md<j> (A.ll) 

J— TV J — TT J — IT J — TT 

for all u. But according to Lemma 5.3, this is only possible if, up to a common rotation, g = g* and h = h±. 
Now repeat the argument for each of the coordinates. A similar argument applies to (2.6). ■ 

Proof of Lemma A.l Given a unit vector u, let S + denote the set of points x in S 1 such that 

(x-u)(u-u ) >0 , (A.12) 

where as before, uo is the "north pole" (0, 1). Let S- be the complement of S + . That is, the line fixed by 
R u slices S* 1 in two, and S + is the half containing the north pole. Then since R u is a measure preserving 
transformation, 

' W g(e)K(\6-cf>\)h(cj>)d6d<j> = 

g(6)K(\0-cj>\)h(<j>)d6d<j>+ [ f g{6)K{\6 - 0|)/i(0)d0d0+ 

S+ JS + JS- JS- 



f f g(0)K(\6-<j>\)h(<j>)d6d(j>+ [ j g(0)K(\6 - 0|)ft(0)d0d<£ = 
f f g(6)K(\e-<j>\)h( ( j))ded ( j)+ f f g(R u 6)K(\6 - <p\)h{R u cp)&6dcp 
f 9 (R u 0)K(\R u 6-cP\)h(<}>)d9d<f>+ f f g{6)K(\0 - i^DM^t 

S+ J S+ J S.\- J S+ 



(A.13) 



The desired inequality is then a consequence of the following inequality for paris of real numbers: Let 
a\ and a 2 and bi and hi be any for positive real numbers. Rearrange a\ and a 2 to decrease, and b\ and 
b 2 to increase; i.e., let Let a\ — max{ai,a 2 } and let a 2 = min{ai,a 2 }, and let 6£ — min{&i,6 2 } and let 
b\ = max{6i,6 2 }. Then 

a\b\ + a* 2 b* 2 < a\b\ + a 2 b 2 , (A. 14) 

and there is equality if and only if a\ = a* and &i = b\ or a\ — a 2 ,b\~b 2 . 
We now apply this with 

a 1 =g(6) a 2 = g{R u 0) h = h{<j>) and b 2 = h(R u (f>) . (A.15) 

Then 

a\ = R+g(6) a* = R+g(R u 0) b\ = R^h(<f>) and b* = R~h(R u cl>) . (A.16) 
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Since 

K(\e - R u <p\) = K{\R u e - <j>\) < K(\e - <f>\) (A.17) 

almost everywhere, we have that 

g(6)K(\e - cf>\)h(<f>) + 9(RuO)K(\6 - <j>\)h{R u cj>)+ 
g(9)K(\R u 9 - cj>\)h{cj>) + g{R u 6)K{\R u 6 - 4>\)RZh{R u <P) > 
R+g(8)K(\8 - 4>\)KK4>) + R+g(R u e)K(\8 - t\)R-h(R u <j>)+ 
R+g(8)K(\R u 8 - <f>\)R^h(<f>) + R+g{R u 8)K{\R u 8 - cj>\)R-h{R u ct>) 



(A.18) 



for almost every 8 and <f> in <S + , with equality if and only if 

g(Ru8) < g(8) and h{R u cf>) > g{<f>) (A. 19) 

or 

g{Ru8)>g(9) and h{R u 4>) < g{<j>) (A.20) 

for almost every 8 and cj> in S+. Now unless g is constant, we can find a 8 and u so that either g(R u 9) < g{9) 
or g(R u 8) > g{9). Suppose it is the first case. Then (A. 19) holds, and for almost every <f), we must have 
h(R u <fi) > g(4>). Making a similar argument for h, we see that one of (A. 19) or (A.20) must hold for almost 
every 9 and <f>. The only way that this can happen is if g and h are symmetric monotone. 
The key here is the following pointwise inequality: If a > b and c > d, then 

F(a, d) + F{b, c) < F(a, c) + F{b, d) . (A.21) 

To see this, let k = c — d and h = a — b. Then 

(F(a,d) + F(b,c))-(F(a,c) + F(b,d)) = -hk [ f 9 f (ft + sh, d + tk)dsdt . (A.22) 

Jo Jo opidp 2 

The second inequality of Lemma A.l follows directly from this. ■ 

Remark As can be seen from the proof, the decrease upon rearrangement can be estimated quantitatively if 
we strengthen (1.2) so that the positive lower bound is uniform. This could be used to relax the conditions 
on the interactions. 



Appendix B. 

In this appendix we prove the following theorem which extends Theorem 2.3 to the case where no a priori 
bound on the densities is assumed, or the entropy term does not necessarily prevents vacuum, or both. 

Theorem B.l Assume: 
i) G and D strictly convex functions 
ii) limbec ^p- = oo , 

Hi) G and D satisfy the following doubling condition: there exists a constant K so that for all r, 

G(2r) < KG(r) and D(2r) < KD{r) . (B.l) 
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iv) there is an L > — oo so that for all r > 0, 

G(r) > L and D{r) > L . (B.2) 

Then, under the same assumptions on the interactions of Theorem 2.3, the conclusions of Theorem 2.3 are 
still true. 

Proof: Wc first observe that for any fixed admissible p^\ the functional Q(p2) defined by 

g{ P2 ) = Hp { i\p2) (b.3) 

is strictly convex on the set of densities satisfying the mean density constraint. By Fatou's lemma, it is also 
lower semicontinuous in the L 1 topology. 

Now since we seek minimizers of J 7 , we may assume that p^ is symmetric monotone about the point on 
the torus that is antipodal to the origin, and then wc may restrict our consideration of p 2 <— ► J-(pf\ p 2 ) to 
densities that are symmetric monotone about the origin. 

Now for any fixed constant M, the set of densities p 2 such that p 2 (x) < M for all x, and that p 2 is 
symmetric monotone is strongly compact in i 1 (A) by the Helly selection principle. It therefore follows that 
there is a unique minimizer p 2 of p 2 JF( / 9^ \p 2 ) in this class. By Lemma 2.2, this is also the unique 
minimizer of Q(p 2 ) in the class of densities satisfying: 

(i) J A p 2 dx = n 2 . 

(ii) p 2 is symmetric monotone with its maximum at the point antipodal to the origin. 

(hi) ||p 2 ||oo <M 

We now show that for M large enough, this minimizer ceases to depend on M, so that the condition (iii) 
becomes superfluous. 

Consider a variation p 2 + h of p 2 . Clearly, we must have h < on 

A M - {x \p 2 (x) = M} , (B.4) 

and ft, > on 

A z - {x \p 2 {x) = 0} . 

We must also have 

h{x)dx = (B.5) 



in order to preserve the condition that (pi , p 2 + h) belongs to V(n\,n 2 ). The Euler-Lagrange condition 
then is 

J (F\p 2 )+D\pf ) +p 2 ) + U*pf ) )hdx>0 (B.6) 

Let B denote (A M U A z f . We hrst show that for M sufficiently large, \B\ ^ 0. Indeed, if \B\ = 0, then 
p 2 = M1a m - Since § A p 2 (x)&x = ni|A|, \Az\ then equals n 2 \A\/M and consequently, 

9(h) > ■ ( B - 7 ) 

Since F(M)/M tends to infinity with M, and since Q(p 2 ) < G(n 2 /\A\) whenever M > n 2 /\A\ so that 
P2 = i2/|A| is an admissible trial density, we obtain a contradiction. Henceforth take M large enough to 
ensure that \B\ ^ 0. Note that on B that the Euler Lagrange equation 

F\p 2 ) + D'{p ( ? ] +p 2 ) + U* pf ] = C (B.8) 
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holds for some value of C, since on this set, if h is an admissible variations, so is —h. 
We have the Chebyshev estimate 



w < n -p ■ <*»> 

Further increasing M, we may assume that \Az\ < |A|/3. It then follows that 

either \B\ > J^l or \A Z \ > ^ . (B.10) 

Our next goal is to obtain an a-priori bound on C. We will obtain two such bounds: one that is valid, when 
\B\ is not too small, and one that is valid when \Az\ is not too small. By the above, at least one of these 
two must hold. 

To obtain a bound that will be useful if \B\ is not too small, integrate the Euler Lagrange equation over 
B, and obtain 

/ F'(p 2 )dx + f D'{pf+p 2 )dx+ f U*pf ] Ax = C\B\ . (B.ll) 

JB JB J B 

Now, by convexity, for any a, 

F'(a) < F(a + 1) - F(a) < F(a + 1) = F (j^j < \F{2a) + \f{2) . (B.12) 
Making use of the doubling condition, we finally have 

F'(a) < |F( a ) + \F{2) . 
Replacing a by p 2 and integrating over B, 

j B F'{p 2 )dx < | jT F(p 2 )dx + if(2)|A| . (B.13) 

In the same way, we obtain 

J &<j>W+f»)dx<^J D(pW+ p 2 )Ax+ l -D{2)\\\ . (B.14) 

Finally, 



IB 

Combining estimates we have 



/ U*p ( i ) dx < a\A\m . (B.15) 

JB 



C, -^l(f e(n2/|A|) + ^ (F(2) + ^ (2))+ani|A| ' ' (B - 1(i) 



so that under the first alternative above, we have the a-priori bound 



C ~ ]Xf (f G ^/^ + \ + D(2)) + ani|A| ) . (B.:i 7 1 
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Next, we obtain an estimate that will be useful is \ Az \ is not too small. First, if \Az\ ^ 0, we can consider 
a variation of the following type: Let h satisfying (B.5) with h + supported by Az, and with h- supported 
in B. Then if we let a = j A h+dx, we have from (B.6) that 

F'(0) + i J (d'(^ 0) ) + U * p { ° ] ) h+dx > C . (B.18) 
By this variational inequality, 

C < F'(0) + [ (d'(p^) + U * p<°>) dx . (B.19) 

\ A Z\ JAz v 7 

The same convexity arguments employed above now yield 

C < F'(0) + j±- ^|c7(n 2 /|A|) + l -D{2) + am|A|) , (B.20) 

which, under the second alternative above, becomes 

C < F'(0) + ± (^Q(n 2 /\A\) + l -D{2) + am|A|) . (B.21) 

Now notice that in the case F(t) = tlogt, F'(0) = —oo, so in this case, \Az\ > is precluded. We therefore 
have the a-priori estimate 

C < max{F'(0), 0} + A (^^(n 2 /|A|) + \ (F(2) + D{2)) + ani \A\) ■ (B.22) 

To apply his, we suppose \Am\ > 0, and consider a variation of the following type: Let h satisfy (B.5) 
with h-, the negative part of h, supported by Am, and with h+ supported in B. Then if we let a — J A h + dx, 
we have from (B.6) that 

-J (f 1 (p 2 ) + D 1 (p^ +p 2 ) + U*p ( °^h^dx + Ca>0 (B.23) 



This means that 



\ J (f'{p 2 ) + D'(p 2 ) + U * /4 0) ) h-dx < C , (B.24) 



I A 

which in turn, since D' is an increasing function, means that 

F'(M) < C . (B.25) 

Let G be the inverse function to F'\ i.e., 

G(t) = mi{r > | F'(r) > t } . (B.26) 

Then, we have the a-priori bound 

M < G(C) , 

which, when combined with our a-priori bound on C implies that for sufficiently large M, \Am\ = 0. 

29 



Therefore, our minimizer p 2 satisfies 



p 2 < G (max{F'(0), 0} + ± G(n 2 /\ A|) + \ (F(2) + D(2)) + am |A|) ) . (B.27) 

Moreover, we see that wherever p 2 ^ 0, it satisfies the Euler-Lagrange equation. 
Now fix this minimizer p 2 and consider the functional 

Qi(pi) = T{ Pu p 2 ) . (B.28) 

Exactly the same argument shows that this has a minimizer pi in the class of densities that are symmetric 
monotone with a maximum at the origin, and with J A pidx — |A|m. Moreover, we obtain in this way the 
a-priori bound 

px < G (max{F'(0), 0} + JL a(m/|A|) + ± (F(2) + D{2)) + an 2 \A\^j ) . (B.29) 

Now consider a minimizing sequence {Pi°\ pj^) in T>(ni,n 2 ) for .T 7 . By the argument above, we can 
replace each pair (p^\p^) by another pair (p^ , p^ ) in D(ni,n 2 ) so that p^ and p^ are symmetric 
monotone about the origin and its antipodal point respectively, and so that p[ and p 2 satisfy the a-priori 
L°° bounds above, and last but not least, so that 

HA k \p2 k) )<Hp ( "\p { 2 k) )- (B.30) 
Now by the Helly selection principle again, we have that for a subsequence, 

pi = lim pf n) (B.31) 

n— >oo 

and 



p 2 = lim pi fen) (B.32) 

exist almost everywhere. By the a-priori L°° bounds, and the dominated convergence theorem, these limits 
also hold in L 1 , and hence (pi,p 2 ) £ T>(m,n 2 ). Moreover, by the lower semicontinuity of T discussed at the 
outset, 

T(p 1 ,p 2 )<limMT(pf"\pf- ) ) = inf ^(pi,p 2 ). (B.33) 

rwoc (pi,P2)eX> (711,712) 

This proves the existence of our minimizers. 

We now drop the tilde, and examine their properties. By the argument above, we clearly have 

Pl (x) =G(G\-U* p 2 (x) - D'( Pl (x) + p 2 {x))) (B.34) 

on the set where pi(x) ^ 0, and likewise 

P2 (x) =G(C 2 -U* prix) - D'( Pl (x) + p 2 {x))) (B.35) 

on the set where p 2 (x) ^ 0, and we have the asserted L°° bounds. Finally, since p 2 and pi are monotone 
symmetric, these L°° bounds imply L 1 bounds for the gradients of pi and p 2 . ■ 
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